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We extend the semiclassical model of the collective atomic recoil laser (CARL) to include the 
quantum mechanical description of the center-of-mass motion of the atoms in a Bose-Einstein con- 
densate (BEC). We show that when the average atomic momentum is less than the recoil momentum 
hq, the CARL equations reduce to the Maxwell-Bloch equations for two momentum levels. In the 
conservative regime (no radiation losses), the quantum model depends on a single collective pa- 
rameter, p, that can be interpreted as the average number of photons scattered per atom in the 
classical limit. When p ^> 1, the semiclassical CARL regime is recovered, with many momentum 
■ levels populated at saturation. On the contrary, when p < 1, the average momentum oscillates 

' between zero and hq, and a periodic train of 2ir hyperbolic secant pulses is emitted. In the dissipa- 

tive regime (large radiation losses) and in a suitable quantum limit, a sequential superfluorescence 
scattering occurs, in which after each process atoms emit a tt hyperbolic secant pulse and populate 
a lower momentum state. These results describe the regular arrangement of the momentum pattern 
observed in recent experiments of superradiant Rayleigh scattering from a BEC. 
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Recently, we have shown 
in terms of the CARL mec 



With the realization of Bose Einstein condensation (BEC) in dilute alkali gases H], it is now possible to study the 
coherent interaction between light and an ensemble of atoms prepared in a single quantum state. For example, Bragg 
diffraction O of a Bose-Einstein condensate by a moving optical standing wave can be used to diffract any fraction 
of a BEC into a selectable momentum state, realizing an atomic beam splitter. Among the multitude of experiments 
, studying the behavior of a BEC under the action of external laser beams, only a small number have been devoted to 
' the active role caused by the atoms in the condensate on the radiation, fn particular, collective light scattering and 
t-H | matter-wave amplification caused by coherent center-of mass motion of atoms in a condensate illuminated by a far 
off- resonant laser were recently observed These experiments have been interpreted in Ref. || as superradiant 

Rayleigh scattering, and successively investigated in Ref. || and using a quantum theory based on a quantum 
multi-mode extention of the collective atomic recoil laser (CARL) Hamiltonian model originally derived by Bonifacio 
et al. P-pd|. fn particular, the original semiclassical CARL model was extended in Ref. (I^JlI to include a quantum 
mechanical description of the center-of-mass motion of the atoms in the condensate. Whereas the analysis of ref. 
p2| , p^[ is limited to the study of the onset of the collective instability starting from quantum fluctuations, some 
Q \ nonlinear effects due to momentum population depletion were discussed in ref. |6[ and . 

14[ | that the superradiant Rayleigh scattering from a BEC can be satisfactorily interpreted 
ranism using a semiclassical model in the 'mean-field' approximation fl5|j , in which the 
rapid escape of the radiation from the condensate is modelled by a decay of the field amplitude at the rate n c — c/2L, 
where L is the sample length. The main drawback of the semiclassical model is that, as it considers the center-of-mass 
motion of the atoms as classical, it cannot describe the discreteness of the recoil velocity, as has been observed in the 
experiment of Ref. j| . 

The aim of this work is to extend the semiclassical CARL model to include the quantum mechanical description 
of the center-of-mass motion of a sample of cold atoms. The quantum model that we obtain is equivalent to that 
derived by Moore and Meystre |Q using second quantization techniques. However, whereas the work of Ref. g| 
is focussed on the linear regime and on the start-up of the instability, we study the full nonlinear regime and the 
quantum and classical limits of the model. Our basic result is that the atomic motion is quantized when the average 
recoil momentum is comparable to hq (where q = k-z — fei is the difference between the incident and the scattered 
wave vectors), i.e. the recoil momentum gained by the atom trading a photon via absorbtion and stimulated emission 
between the incident and scattered waves. In this limit, the quantum CARL equations reduce to the Maxwell-Bloch 
equations for two momentum levels |l6f| . In the 'conservative' (or 'hamiltonian') regime, in which the radiation losses 
are negligible, this occurs for p < 1, where the CARL parameter p represents the average number of photons scattered 
per atom in the classical limit. In the superradiant regime, for k > f (where k = n c /uj r p, k c is the radiation loss, w r p is 
the collective recoil bandwidth, ui r — h\q\ 2 /2M is the recoil frequency and M is the atomic mass), the atomic motion 
becomes quantized for p < \/2k. In this limit, we demonstrate that a sequential superfluorescence (SF) scattering 
occurs, in which, during each process, the atoms emit a 7r hyperbolic secant pulse and populate a lower momentum 
level, as it has been observed in the MIT experiment || . 
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Our starting point is the classical model of equations for N two-level atoms exposed to an off-resonant pump laser, 
whose electric field Eq = eSg cos(fc2 • x — LU2t) is polarized along e, propagates along the direction of ki and has 
a frequency 102 — ck 2 with a detuning from the atomic resonance, A20 = u>2 — <-^o, much larger than the natural 
linewidth of the atomic transition, 7. We assume the presence of a scattered field ('probe beam') with frequency 

lo\ == L02 A21, wave-number k\ making an ang le (j) with k 2 and electric field E = (e/2)[£ (tje^ 1 ^-^ + c.c] with the 
same polarization of the pump field. In the absence of an injected probe field, the emission starts from fluctuations and 
the propagation direction of the scattered field is determined either by the geometry of the condensate (as in the case 
of the MIT experiment |^] , where the condensate has a cigar shape) or by the presence of an optical resonator tuned 
on a selected longitudinal mode. By adiabatically eliminating the internal atomic degree of freedom, the following 
semiclassical CARL equations has been derived p-^0|: 



(1) 
(2) 



or L 

A A N 

3=1 

where r = poj r t is the interaction time in units of u> r p, 9j = (ki — k 2 ) ■ xj = qzj and pj — qv Z j / pu r (where q = \q\) are 
the dimensionless position and velocity of the j'-th atom along the axis z, A — —i(eo/n s hu>p) 1 / 2 £(T)e tST , 5 = A 2 i/cJ r p 
and p = (Qo/2A2o) 2 / 3 (up 2 n s /heoLUr) 1 / 3 is the collective CARL parameter, flo = p£o/h is the Rabi frequency of the 
pump, n s = N/V is the average atomic density of the sample (containing N atoms in a volume V) and p, is the dipole 
matrix element. We assume u>2 ~ — ck, so that q m 2k sin 2 (^/2). Eqs.(|l])-(^) are formally equivalent to those of 
the free electron laser model ]ri| . In order to quantize both the radiation field and the center-of-mass motion of the 
atoms, we consider 9j, pj = (p/2)pj — Mv Z j/hq and a = (N p/2) 1 / 2 A as quantum operators satisfying the canonical 
commutation relations [0j,Pj'] — iSjj' and [a, at] = 1. With these definitions, Eqs.(|])-(||) are the Heisenberg equations 
of motion associated with the Hamiltonian: 



1 N 



N IN \ N 

P j=l \j=l J 3=1 



where g = y/ p/2N. We note that [H, Q] — 0, where Q = a^a + J2j=iPj 1S ^ ne total momentum in units of hq. 
In order to obtain a simplified description of a BEC as a system of N noninteracting atoms in the ground state, 
we use the Schrodinger picture for the atoms (instead of the usual Heisenberg picture Jl8|), i.e. \tp(0i, . . . ,6n)) = 
\i(>(9i)) . . . \ip(6ff)), where \ip(0j)) obeys the single-particle Schrodinger equation, i(d / f dT)\ip(9j)) = Hj(6j,p : j)\'t/j(0j)). 
In this paper we describe the light field classically. Hence, considering the field operator a as a c- number, eq.(H) yields: 

, N 

£ =l 5a + gJ2(m)\e-^\m))- (5) 

Let now expand the single-atom wavefunction on the momentum basis, \ip(0j)) = Y]„ Cj(n)\n}j, where pj \n)j = n\n)j, 
n = —00, . . . 00 and Cj(n) is the probability amplitude of the j-th atom having momentum —nhq. Introducing the 
collective density matrix: 

1 N 

s m , n =^Y. c ^ m T c ^ n y {m ~ n)s ^ ( 6 ) 

a straightforward calculation yields, from Eqs.(^) and (||), the following closed set of equations: 

, m '" = i(m - n)5 m ,„5 m , n + ^ [A (S m +i, n ~ S m ,n-i) + A* (S m , n +i - S m -i,n)] (7) 

(XT 2 

dA x 

-j- = S n>n+ x - kA, (8) 

n= — 00 
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where 5 m . n = 5+(m + n)/p and we have redefined the field as A = yj2/ pNae~ l5T . We have also introduced a damping 
term —kA in the field equation, where n = n c /uj r p, k c = c/2L and L is the sample length along the probe propagation, 
which provides an approximated model describing the escape of photons from the atomic medium. In the presence of 
a ring cavity of length L cav and reflectivity R, n c = — (c/L cav )lni?, as shown in the usual 'mean-field' approximation 
jn]]. Eqs.(0) and (||) are completely equivalent to the CARL equations (|l|)-(j^) and determine the temporal evolution 
of the density matrix elements for the momentum levels. In particular, p n — S„ tn is the probability of finding the atom 
in momentum level |n), (p) = Y] n nS n ,n is the average momentum and ^ n SVi,n+i is the bunching parameter. Eqs.(0) 
and (^) are identical to these derived by Moore and coworkers [fl2| second quantizing the single-particle Hamiltonian 
Hj and introducing bosonic creation and annihilation operators of a given center-of-mass momentum. For a constant 
field A, Eq.(^|) describes a Bragg scattering process, in which m — n photons are absorbed from the pump and scattered 
into the probe, changing the initial and final momentum states of the atom from m to n. Conservation of energy and 
momentum require that during this process ui\ — u>2 = (m + n)ui r , i.e. 5 m ^ n — 0. Eqs.(^) and (^) conserve the norm, 
i- e - Em Srn,m = 1, and, when k = 0, also the total momentum Q = (p/2)\A\ 2 + (p). Figure [l]a shows \A\ 2 vs. r, for 
k = 0, S = and A(0) = 1CP 4 , comparing the semiclassical solution with the quantum solution in the classical limit, 
p ^> 1: the dashed line is the numerical solution of Eqs.(0)-(^), for a classical system of N = 200 cold atoms, with 
initial momentum Pj(0) — (where j = 1, . . . , N) and phase 0j(O) uniformly distributed over 2tt 7 i.e. unbunched; 
the continuous line is the numerical solution of Eqs.(0) and (0) for p — 10 and a quantum system of atoms initially 
in the ground state n = 0, i.e. with S n>m = S n oS m o. Figure Ha shows that the quantum system behaves, with good 
approximation, classically. Because from Figjlja the maximum dimensionless intensity is \A\ 2 ps 1.4, the constant of 
motion Q gives (p) rs —0.7/9 and the maximum average number of emitted photons is about (era) ~ N p. Hence, the 
CARL parameter p can be interpreted as the maximum average number of photons emitted per atom (or equivalently, 
as the maximum average momentum recoil, in units of Kq, acquired by the atom) in the classical limit. Figure ^> 
shows the distribution of the population level p n at the first peak of the intensity of figj^a, for r = 12.4. We observe 
that, at saturation, twenty-five momentum levels are occupied, with an induced momentum spread comparable to the 
average momentum. 

Let us now consider the equilibrium state with no probe field, A — 0, and all the atoms in the same momentum 
state n, i.e. with SVi,n = 1 and the other matrix elements zero. This is equivalent to assume the temperature of the 
system equal to zero and all the atoms moving with the same velocity —nhq, without spread. This equilibrium state 
is unstable for certain values of the detuning. In fact, by linearizing Eqs.(^) and (^) around the equilibrium state, 
the only matrix elements giving linear contributions are S n —i tn and S^n+ii showing that in the linear regime the 
only transitions allowed from the state n are these towards the levels n — 1 and n + 1. Introducing the new variables 
B n = Sn.n+i + Sn-i.n and P n = S n ,n+i ~ S n -i,n> Eqs.(0) and (||) reduce to the linearized equations: 

(9) 

P A (10) 

(11) 

where S n = S + 2n/p. Seeking solutions proportional to e l< - A_<5 ™- ) ' r , we obtain the following cubic dispersion relation: 

(A — S n — in) (A 2 — 1/p 2 ) + 1 = 0. (12) 

In the exponential regime, when the unstable (complex) root A dominates, B(t) ~ e J ( A -< 5 ») T and, from Eq.(^|), 
P n = —p\B n . The semiclassical limit is recovered for p~^> 1 (when k = 0) or p ^> y/H (when k > 1) and S n » 5, i.e. 
neglecting the shift due to the recoil frequency u> r . In this limit, maximum gain occurs for S = 0, with A = (l-n/3)/2 
when k = or A = —(1 + i)/y2~K when k > 1. Furthermore, |iSVi,n+l| ~ l^n-i.nl, so that the atoms may experience 
both emission and absorbtion. This result can be interpreted in terms of single-photon emission and absorpion by 
an atom with initial momentum —nhq. In fact, energy and momentum conservation impose lo\ — L02 — (2n =p l)tv r 
(i.e. 5 n = ±1/(0 ) when a probe photon is emitted or absorbed, respectively. Because in the semiclassical limit the 
gain bandwidth is Aw ~ co r p 3> oj r when k — (or Aw ~ k c 3> lu t when k > 1) the atom can both emit or 
absorbe a probe photon. On the contrary, in the quantum limit the recoil energy huj r can not be neglected, and 
there is emission without absorbtion if |5 n , n +i| -C |5 n _i in |, i.e. B n w —P n and A s» 1/p. This is true for p < 1 
when k = 0, with the unstable root Aw 1/p + 5' n /2 — l/2y/8 12 — 2p (where S' n = S n — 1/p), and for p < \[2k when 
K > 1, with RcA l/p + (pS' n /2)/(5' 2 + k 2 ) and ImA « —{pn/2)/{8 12 + k 2 ). In both cases, maximum gain occurs for 
8 n = 1/p (i.e. A21 = (1 — 2n)uj r ) within a bandwidth Alu ~ uj r p 3 ' 2 and Aco ~ tu r p 2 / n (respectively for k = and 
k > 1), which arc both less than the frequency difference 2cu r between the emission and absorbtion lines. Hence, in 



dB n 
dr 

dP n 
dr 



2 S rii Pn 



= _ 7 A p B 

tun* 7-1 J-^r, 



dA 

—— — D n — KA, 

dr 
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the quantum limit the optical gain is due exclusively to emission of photons, whereas in the semiclassical limit gain 
results from a positive difference between the average emission and absorbtion rates. When k = 0, the resonant gain 
in the limit p < 1 is Gs — u> r p^p/2 = y / 3787r(^o/2A2o)7v / ^eff> where 7 = p 2 k 3 /3Trhea is the natural decay rate 
of the atomic transition and N e s = (A 2 /A)(c/ r yL)N is the effective atomic number in the volume V = AL, where A 
and L are the cross section and the length of the sample. When k > 1, the resonant SF gain in the limit p < \J~2k~ is 
G SF = w rP 2 /2n = (3/ATr) 1 (n /2A 20 ) 2 (X 2 /A)N. 

The above results show that the combined effect of the probe and pump fields on a collection of cold atoms in 
a pure momentum state n is responsible of a collective instability that leads the atoms to populate the adjacent 
momentum levels n — 1 and n + 1. However, in the quantum limit p < 1 when k = (or p < \p2~k~ when k > 1) 
conservation of energy and momentum of the photon constrains the atoms to populate only the lower momentum 
level n — 1. This holds also in the nonlinear regime, as we have verified solving numerically Eqs.(^) and (||). In the 
quantum limit above, the exact equations reduce to those for only three matrix elements, S n<n , S n _i : „_x and S„_i „, 
with >S„-i. „-i + S n>n = 1. Introducing the new variables S n — S n _x >n and W n — S n>n — S'n-i.n-i) Eqs.(Q) and (g) 
reduce to the well-known Maxwell-Bloch equations |L9| : 

^ = - iS ' n S n + P -AW n (13) 

(XT 2 



dWn 

dr 
dA 

~dr~ 



-p(A*S n + h.c.) (14) 
S n - kA, (15) 



where S' n = 5 + (2n — l)/p. When k = and S' n = 0, if the system starts radiating incoherently by pure quantum- 
mechanical spontaneous emission, the solution of Eqs.(p^)-(p^|) is a periodic train of 2tt hyperbolic secant pulses |2C| ] 
with \A\ 2 = (2/p)sech 2 [ Ak / p/2(r — r„)], where r„ = (2n + l)ln(p/2)/-^/ p/2. Furthermore, the average momentum 
(p) = n + Th 2 [yj p/2(r — T n )] — 1 oscillates between n and n — 1 with period r n . We observe that the maximum 
number of photons emitted is (a^a) pca k = {pN/2)\A\ 2 cak = N, as expected. Figured shows the results of a numerical 
integration of Eqs.(0) and (]§), for k = 0, p = 0.2 and 5 = 5, with the atoms initially in the momentum level n = 
and the field starting from the seed value Aq = 10~ 5 . Figures ||a and b show the intensity \A\ 2 and the average 
momentum (p) vs. r, in agreement with the predictions of the reduced Eqs.(p^)-([T5|). 



In the superradiant regime, k > 1, Eqs.(13)-([L5|) describe a single SF scattering process in which the atoms, 
initially in the momentum state n, 'decay' to the lower level n — 1 emitting a n hyperbolic secant pulse, with intensity 
\A\ 2 = 1/[4(k 2 + <5' 2 )]sech 2 [(r - t d )/t sf ] and average momentum (p) = n - (1/2){1 + Th[(r - t d )/t sf }}, where 



tsf = 2(k 2 +5' 2 )/ pn is the 'superfluorescence time' 15|, td = T5FArcsech(2|5'„(0)|) ~ — TsF^\/2\S n (0)\ is the delay 
time and |5 n (0)| <C 1 is the initial polarization. Figures |^a and b shows \A\ 2 and (p) vs. r calculated solving Eqs.(|7]) 
and (||) numerically with n — 10, p = 2, S = 0.5 and the same initial conditions of fig. ||. We observe a sequential 
SF scattering, in which the atoms, initially in the level n = 0, change their momentum by discrete steps of fiq and 
emit a SF pulse during each scattering process. We observe that for 5 = 1/p the field is resonant only with the first 
transition, from n = to n = —1; for a generic initial state n, resonance occurs when S = (1 — 2n)/p, so that in the 
case of figure ||a the peak intensity of the successive SF pulses is reduced (by the factor l/[n 2 + (2n/p) 2 }) whereas the 
duration and the delay of the pulse are increased. However, the pulse retains the characteristic sech 2 shape and the 
area remains equal to 7r, inducing the atoms to decrease their momentum by a finite value Tiq. We note that, although 
the SF time in the quantum limit (tsf = 2k/ p at resonance) can be considerable longer than the characteristic 
superradiant time obtained in the classical limit, tsr = V2k, the peak intensity of the pulse in the quantum limit is 
always approximately half of the value obtained in the semiclassical limit (see Ref. Jlij for details). 

In conclusion, we have shown that the CARL model describing a system of atoms in their momentum ground state 
(as those obtained in a BEC) and properly extended to include a quantum-mechanical description of the center-of-mass 
motion, allows for a quantum limit in which the average atomic momentum changes in discrete units of the photon 
recoil momentum fiq and reduce to the Maxwell-Bloch equations for two momentum levels. These results demonstrate 
that the regular arrangement of momentum pattern observed in the MIT experiment || can be interpreted as being 
due to sequential superfluorescence scattering. A detailed study of this and other aspects of the MIT experiment will 
be the object of a future extended pubblication. 

We thank G.R.M. Robb for a careful reading of the manuscript and for helpful discussions. 
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FIG. 1. Classical limit of CARL for p S> 1 in the case k = 0. (a): \A\ 2 vs. r as obtained from the classical eqs.(l)-(3) 
(dashed line) and from the quantum eqs.(7) and (8) for p — 10 (solid line); (b): population level p n vs. n at the occurring of 
the first maximum of \A\ 2 , at r = 12.4. The other parameters are 5 = and A(0) = 10~ 4 . 



FIG. 2. Quantum limit of CARL for p < 1 in the case n = 0. (a) \A\ 2 and (b) (p) vs. r, for p = 0.2, 5 = 5, A(0) = 10~ 5 and 
the atoms initially in the state n = 0. We note that (p) = -{2/p)(\A\ 2 - \A(0)\ 2 ). 



FIG. 3. Sequential superfluorescent (SF) regime of CARL, (a) \A\ 2 and (b) (p) vs. r, for p = 2, 5 = 0.5, k = 10, and the 
same initial conditions of fig. 2. 



5 



Quantum effects in the collective light scattering by coherent atomic recoil in a 

Bose-Einstein condensate 

N. Piovella, M. Gatelli and R. Bonifacio 
Dipartimento di Fisica, Universita Degli Studi di Milano, INFN & INFM 
Via Celona 16, Milano 1-20133, Italy 

We extend the semiclassical model of the collective atomic recoil laser (CARL) to include the 
quantum mechanical description of the center-of-mass motion of the atoms in a Bose-Einstein con- 
densate (BEC). We show that when the average atomic momentum is less than the recoil momentum 
hq, the CARL equations reduce to the Maxwell-Bloch equations for two momentum levels. In the 
conservative regime (no radiation losses), the quantum model depends on a single collective pa- 
rameter, p, that can be interpreted as the average number of photons scattered per atom in the 
classical limit. When p ^> 1, the semiclassical CARL regime is recovered, with many momentum 
levels populated at saturation. On the contrary, when p < 1, the average momentum oscillates 
between zero and hq, and a periodic train of 2ir hyperbolic secant pulses is emitted. In the dissipa- 
tive regime (large radiation losses) and in a suitable quantum limit, a sequential superfluorescence 
scattering occurs, in which after each process atoms emit a ir hyperbolic secant pulse and populate 
a lower momentum state. These results describe the regular arrangement of the momentum pattern 
observed in recent experiments of superradiant Rayleigh scattering from a BEC. 

PACS numbers: 42. 50. Fx, 42.50.Vk, 03.75.Fi 

With the realization of Bose Einstein condensation (BEC) in dilute alkali gases [1], it is now possible to study the 
coherent interaction between light and an ensemble of atoms prepared in a single quantum state. For example, Bragg 
diffraction [2] of a Bose-Einstein condensate by a moving optical standing wave can be used to diffract any fraction 
of a BEC into a selectable momentum state, realizing an atomic beam splitter. Among the multitude of experiments 
studying the behavior of a BEC under the action of external laser beams, only a small number have been devoted to 
the active role caused by the atoms in the condensate on the radiation. In particular, collective light scattering and 
matter-wave amplification caused by coherent center-of mass motion of atoms in a condensate illuminated by a far 
off-resonant laser were recently observed [3-5]. These experiments have been interpreted in Ref. [3] as superradiant 
Rayleigh scattering, and successively investigated in Ref. [6] and [7] using a quantum theory based on a quantum 
multi-mode extention of the collective atomic recoil laser (CARL) Hamiltonian model originally derived by Bonifacio 
et al. [8-11]. In particular, the original semiclassical CARL model was extended in Ref. [12,13] to include a quantum 
mechanical description of the center-of-mass motion of the atoms in the condensate. Whereas the analysis of ref. 
[12,13] is limited to the study of the onset of the collective instability starting from quantum fluctuations, some 
nonlinear effects due to momentum population depletion were discussed in ref. [6] and [7]. 

Recently, we have shown [14] that the superradiant Rayleigh scattering from a BEC can be satisfactorily interpreted 
in terms of the CARL mechanism using a semiclassical model in the 'mean-field' approximation [15], in which the 
rapid escape of the radiation from the condensate is modelled by a decay of the field amplitude at the rate k c = c/2L, 
where L is the sample length. The main drawback of the semiclassical model is that, as it considers the center-of-mass 
motion of the atoms as classical, it cannot describe the discreteness of the recoil velocity, as has been observed in the 
experiment of Ref . [3]. 

The aim of this work is to extend the semiclassical CARL model to include the quantum mechanical description 
of the center-of-mass motion of a sample of cold atoms. The quantum model that we obtain is equivalent to that 
derived by Moore and Meystre [12] using second quantization techniques. However, whereas the work of Ref. [12] 
is focussed on the linear regime and on the start-up of the instability, we study the full nonlinear regime and the 
quantum and classical limits of the model. Our basic result is that the atomic motion is quantized when the average 
recoil momentum is comparable to hq (where q = k'l — k\ is the difference between the incident and the scattered 
wave vectors), i.e. the recoil momentum gained by the atom trading a photon via absorbtion and stimulated emission 
between the incident and scattered waves. In this limit, the quantum CARL equations reduce to the Maxwell-Bloch 
equations for two momentum levels [16]. In the 'conservative' (or 'hamiltonian') regime, in which the radiation losses 
are negligible, this occurs for p < 1, where the CARL parameter p represents the average number of photons scattered 
per atom in the classical limit. In the superradiant regime, for k > 1 (where k = K c /ui r p, k c is the radiation loss, u> r p is 
the collective recoil bandwidth, oj r = h\q\ 2 /2M is the recoil frequency and M is the atomic mass), the atomic motion 
becomes quantized for p < \/2k. In this limit, we demonstrate that a sequential superfluorescence (SF) scattering 
occurs, in which, during each process, the atoms emit a 7r hyperbolic secant pulse and populate a lower momentum 
level, as it has been observed in the MIT experiment [3]. 
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Our starting point is the classical model of equations for N two-level atoms exposed to an off-resonant pump laser, 
whose electric field Eo = e£o cos(&2 • % — <+>2t) is polarized along e, propagates along the direction of k 2 and has 
a frequency u>2 = ck 2 with a detuning from the atomic resonance, A20 = lo 2 — ^0, much larger than the natural 
linewidth of the atomic transition, 7. We assume the presence of a scattered field ('probe beam') with frequency 
lu\ = uj'2 — A21, wavenumber k\ making an angle <j> with £2 and electric field E = (e/2)[£(t)e l ( kl ' x ~ Wlt ) + ex.] with the 
same polarization of the pump field. In the absence of an injected probe field, the emission starts from fluctuations and 
the propagation direction of the scattered field is determined either by the geometry of the condensate (as in the case 
of the MIT experiment [3], where the condensate has a cigar shape) or by the presence of an optical resonator tuned 
on a selected longitudinal mode. By adiabatically eliminating the internal atomic degree of freedom, the following 
semiclassical CARL equations has been derived [8-10]: 

j=i 

where r = pu> r t is the interaction time in units of u> r p, 9j = {k\ — k 2 ) • %j = QZj and fj = Qv z j / pui r (where q = \q\) are 
the dimensionless position and velocity of the j-th atom along the axis z, 

and p = (Slo/2A2o) 2/ ' 3 (w / u 2 n s /7}eowj?) 1/ ' 3 is the collective CARL parameter. Slo = p£o/Ti is the Rabi frequency of the 
pump, n s = N/V is the average atomic density of the sample (containing N atoms in a volume V) and p is the dipole 
matrix element. We assume W2 k, lu\ = ck, so that q m 2k sm 2 ((f)/2). Eqs.(l)-(3) are formally equivalent to those of 
the free electron laser model [17]. In order to quantize both the radiation field and the center-of-mass motion of the 
atoms, we consider 9j, pj = (p/2)p- = Mv z j/hq and a = (N p/2) 1 / 2 A as quantum operators satisfying the canonical 

commutation relations [0j,Pj>] = i6jj> and [a, aJ] = 1. With these definitions, Eqs.(l)-(3) are the Heisenberg equations 
of motion associated with the Hamiltonian: 



jv / N \ N 

P i=i \i=i / i=i 

where g = y/ p/2N . We note that [H , Q] = 0, where Q = a^a + X2j=ift' ^ s ^ e total momentum in units of Tiq. 
In order to obtain a simplified description of a BEC as a system of N noninteracting atoms in the ground state, 
we use the Schrodinger picture for the atoms (instead of the usual Heisenberg picture [18]), i.e. \ip(6i, . . . , 6^)} = 
\tl>(9i)) . . . \tl>(9^)), where \ip(9j)} obeys the single-particle Schrodinger equation, i(d/dT)\ip(6j)) = Hj(9j , Pj)\ip(@j)} ■ 
In this paper we describe the light field classically. Hence, considering the field operator a as a c-number, eq.(3) yields: 

fl N 

£ = i 6 a + 9 J2W0:)\e-^m0:))- (5) 
i=i 

Let now expand the single-atom wavefunction on the momentum basis, \ip(9j)} = ^2 n Cj(n)\n)j, where Pj\n)j = n\n)j, 
n = — cxd, . . . 00 and Cj(n) is the probability amplitude of the j-th atom having momentum —nhq. Introducing the 
collective density matrix: 

1 N 

s m ,n = ^j2 c i( m y c i( n v (m ~ n)ST > ( 6 ) 

i=i 

a straightforward calculation yields, from Eqs.(4) and (5), the following closed set of equations: 

—j—^ = i(/77 ^)*^m,n S m>n -\- — [A (Sm + i n S m>n — i ) + A (iSVr^n + 1 S m — i >n )] (7) 

dA °° 

^ S n ,n + 1 ~ KA, (8) 



dT 



n = — 00 
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where h m n = 8 + (m + n)/ p and we have redefined the field as A = \J2/ pNae~ lSr . We have also introduced a damping 
term — kA in the field equation, where k = K c /ui r p, k c = c/2L and L is the sample length along the probe propagation, 
which provides an approximated model describing the escape of photons from the atomic medium. In the presence of 
a ring cavity of length L cav and reflectivity R, k c = — (c/_L cav )lni?, as shown in the usual 'mean-field' approximation 
[11]. Eqs.(7) and (8) are completely equivalent to the CARL equations (l)-(3) and determine the temporal evolution 
of the density matrix elements for the momentum levels. In particular, p n = S nt n is the probability of finding the atom 
in momentum level \n), (p) = ^2 n nS nt n is the average momentum and ^2 n S n>n +i is the bunching parameter. Eqs.(7) 
and (8) are identical to these derived by Moore and coworkers [12] second quantizing the single-particle Hamiltonian 
Hj and introducing bosonic creation and annihilation operators of a given center-of-mass momentum. For a constant 
field A, Eq.(7) describes a Bragg scattering process, in which m — n photons are absorbed from the pump and scattered 
into the probe, changing the initial and final momentum states of the atom from m to n. Conservation of energy and 
momentum require that during this process u>\ — u>2 = (m + n)u> r , i.e. 8 m n = 0. Eqs.(7) and (8) conserve the norm, 
i.e. Sm^m = 1, and, when k = 0, also the total momentum Q = (/o/2)|j4| 2 + (p). Figure la shows |j4| 2 vs. t, for 
k = 0, 6 = and A(0) = 10 -4 , comparing the semiclassical solution with the quantum solution in the classical limit, 
p 1: the dashed line is the numerical solution of Eqs.(l)-(3), for a classical system of N = 200 cold atoms, with 
initial momentum Pj(0) = (where j = 1,...,N) and phase 0j(0) uniformly distributed over 2tt, i.e. unbunched; 
the continuous line is the numerical solution of Eqs.(7) and (8) for p = 10 and a quantum system of atoms initially 
in the ground state n = 0, i.e. with S nt m = 6 n o8 m o- Figure la shows that the quantum system behaves, with good 
approximation, classically. Because from Fig. la the maximum dimensionless intensity is |j4| 2 ps 1.4, the constant of 
motion Q gives (p) ps —0.7/9 and the maximum average number of emitted photons is about (a^a) ~ Np. Hence, the 
CARL parameter p can be interpreted as the maximum average number of photons emitted per atom (or equivalently, 
as the maximum average momentum recoil, in units of Tiq, acquired by the atom) in the classical limit. Figure lb 
shows the distribution of the population level p n at the first peak of the intensity of fig. la, for r = 12.4. We observe 
that, at saturation, twenty-five momentum levels are occupied, with an induced momentum spread comparable to the 
average momentum. 

Let us now consider the equilibrium state with no probe field, A = 0, and all the atoms in the same momentum 
state n, i.e. with S nt n = 1 and the other matrix elements zero. This is equivalent to assume the temperature of the 
system equal to zero and all the atoms moving with the same velocity —nkq, without spread. This equilibrium state 
is unstable for certain values of the detuning. In fact, by linearizing Eqs.(7) and (8) around the equilibrium state, 
the only matrix elements giving linear contributions are S n — \^ n and S n ^ n -\.\ , showing that m the linear regime the 
only transitions allowed from the state n are these towards the levels n — 1 and n + 1. Introducing the new variables 
B n = S n>n +i + S n -\ t n and P n = S n>n +i — S n -i >n , Eqs.(7) and (8) reduce to the linearized equations: 

= _,•*„ B„-ip„ ( 9 ) 

dr p 
dP i 

-T^=-i6 n Pn--B n -pA (10) 

dr p 
dA 

-y = B n — kA, (11) 
dr 

where 6 n = 6 + In/ p. Seeking solutions proportional to e l ( x ~ Sn ' >T , we obtain the following cubic dispersion relation: 

(A - 6 n - in) (A 2 - l/p 2 ) + 1 = 0. (12) 

In the exponential regime, when the unstable (complex) root A dominates, B(t) ~ e t ( x - s n) T and, from Eq.(9), 
P n = —p\B n . The semiclassical limit is recovered for p 1 (when k = 0) or p s/k (when k > 1) and 6 n ps 6, i.e. 
neglecting the shift due to the recoil frequency u> r . In this limit, maximum gain occurs for 6 = 0, with A = (l-iV3)/2 
when k = or A = —(1 + i)/\/2~K when n > 1. Furthermore, IS^n+iJ ~ l^n-i^l, so that the atoms may experience 
both emission and absorbtion. This result can be interpreted in terms of single-photon emission and absorpion by 
an atom with initial momentum —nkq. In fact, energy and momentum conservation impose u>\ — u>2 = (2n =F l)w r 
(i.e. 6 n = H/p) when a probe photon is emitted or absorbed, respectively. Because in the semiclassical limit the 
gain bandwidth is Aui ~ u> r p ^> u> r when k = (or Aui ~ k c u> r when k > 1) the atom can both emit or 
absorbe a probe photon. On the contrary, in the quantum limit the recoil energy Tiu> r can not be neglected, and 
there is emission without absorbtion if IS^n+iJ <€i \S n -i,n\, i- e - B n ps — P n and A ps l/p. This is true for p < 1 
when k = 0, with the unstable root A ps l/p + S' n /2 — \/2\JbQ — 2p (where 6' n = 6 n — l/p), and for p < \/2k when 
k > 1, with ReA p» l/p+ (p6' n /2)/(8' 2 + k 2 ) and ImA p» —(pn/2)/(8' 2 + k 2 ). In both cases, maximum gain occurs for 
6 n = l/p (i.e. A21 = (1 — 2n)uj r ) within a bandwidth Aw ~ Lo r p 3 l 2 and Aw ~ w r p 2 /k (respectively for k = and 
k > 1), which are both less than the frequency difference 2u> r between the emission and absorbtion lines. Hence, in 
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the quantum limit the optical gain is due exclusively to emission of photons, whereas in the semiclassical limit gain 
results from a positive difference between the average emission and absorbtion rates. When k = 0, the resonant gain 
in the limit p < 1 is Gs = ^> r P\/ p/2 = ^/3/87r(£lo /2A2o)tV N e ff, where 7 = p 2 k 3 /Sirkeo is the natural decay rate 
of the atomic transition and N e ff = (A 2 /A)(c/jL)N is the effective atomic number in the volume V = AL, where A 
and L are the cross section and the length of the sample. When n > 1, the resonant SF gain in the limit p < \/2k is 
G SF = u r p 2 /2n = (3/47r) 7 (fio/2A 20 ) 2 (A 2 M)#. 

The above results show that the combined effect of the probe and pump fields on a collection of cold atoms in 
a pure momentum state n is responsible of a collective instability that leads the atoms to populate the adjacent 
momentum levels n — 1 and n + 1. However, in the quantum limit p < 1 when n = (or p < \/2k when n > 1) 
conservation of energy and momentum of the photon constrains the atoms to populate only the lower momentum 
level n — 1. This holds also in the nonlinear regime, as we have verified solving numerically Eqs.(7) and (8). In the 
quantum limit above, the exact equations reduce to those for only three matrix elements, S nt n, S n -\ t n-i and S n -\ t n, 
with S n -\ t n-i + S n>n = 1. Introducing the new variables S n — S n — \ n and W n — S n n — S n — \ n — ii 

Eqs.(7) and (8) 

reduce to the well-known Maxwell-Bloch equations [19]: 



dS, 

"dT 

dW, 



dT -KS n + P -AW n (13) 



-p(A*S n +h.c.) (14) 
(It 

d A 

-T = S n - kA, (15) 
(It 

where 6' n = 6 + (2n — l)/p. When k = and 6' n = 0, if the system starts radiating incoherently by pure quantum- 
mechanical spontaneous emission, the solution of Eqs.(13)-(15) is a periodic train of 2tt hyperbolic secant pulses [20] 
with 1 2 = (2/ / o)sech 2 [ v / p72(r - r n )], where r n = (2n + l)ln( / o/2)/ ' \fpj2- Furthermore, the average momentum 
(j>) = n + Th 2 [^/ p/2(t — r n )] — 1 oscillates between n and n — 1 with period r n . We observe that the maximum 
number of photons emitted is (a) ^} pe ak = |yl|p eak = N , as expected. Figure 2 shows the results of a numerical 

integration of Eqs.(7) and (8), for k = 0, p = 0.2 and 6 = 5, with the atoms initially in the momentum level n = 
and the field starting from the seed value A§ = 10 -5 . Figures 2a and b show the intensity |j4| 2 and the average 
momentum (p) vs. r, in agreement with the predictions of the reduced Eqs.(13)-(15). 

In the superradiant regime, k > 1, Eqs.(13)-(15) describe a single SF scattering process in which the atoms, 
initially in the momentum state n, 'decay' to the lower level n — 1 emitting a 7r hyperbolic secant pulse, with intensity 
|j4| 2 = 1/[4(k 2 + i5^ 2 )]sech 2 [(r — td)/tsf] an d average momentum (p) = n — (1/2){1 + Th[(r — td)/tsf]}, where 
tsf = 2(k 2 + 6'n)/pK is the 'superfluorescence time' [15], td = TsFArcsech(2|5' n (0)|) m — tsf ln^/2|5' n (0)| is the delay 
time and IS^O)! <C 1 is the initial polarization. Figures 3a and b shows |j4| 2 and (p) vs. r calculated solving Eqs.(7) 
and (8) numerically with k = 10, p = 2, 6 = 0.5 and the same initial conditions of fig. 2. We observe a sequential 
SF scattering, in which the atoms, initially in the level n = 0, change their momentum by discrete steps of Tiq and 
emit a SF pulse during each scattering process. We observe that for 6 = l/p the field is resonant only with the first 
transition, from n = to n = —1; for a generic initial state n, resonance occurs when 6 = (1 — 2n)/p, so that in the 
case of figure 3a the peak intensity of the successive SF pulses is reduced (by the factor 1/[k 2 + (2n/ p) 2 ]) whereas the 
duration and the delay of the pulse are increased. However, the pulse retains the characteristic sech 2 shape and the 
area remains equal to 7r, inducing the atoms to decrease their momentum by a finite value Tiq. We note that, although 
the SF time in the quantum limit (tsf = 2k/ p at resonance) can be considerable longer than the characteristic 
superradiant time obtained in the classical limit, tsr = \/2k, the peak intensity of the pulse in the quantum limit is 
always approximately half of the value obtained in the semiclassical limit (see Ref. [14] for details). 

In conclusion, we have shown that the CARL model describing a system of atoms in their momentum ground state 
(as those obtained in a BEC) and properly extended to include a quantum-mechanical description of the center-of-mass 
motion, allows for a quantum limit in which the average atomic momentum changes in discrete units of the photon 
recoil momentum Tiq and reduce to the Maxwell-Bloch equations for two momentum levels. These results demonstrate 
that the regular arrangement of momentum pattern observed in the MIT experiment [3] can be interpreted as being 
due to sequential superfluorescence scattering. A detailed study of this and other aspects of the MIT experiment will 
be the object of a future extended pubblication. 

We thank G.R.M. Robb for a careful reading of the manuscript and for helpful discussions. 
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FIG. 1. Classical limit of CARL for p ^> 1 in the case k = 0. (a): |y!| 2 vs. r as obtained from the classical eqs.(l)-(3) 
(dashed line) and from the quantum eqs.(7) and (8) for p = 10 (solid line); (b): population level p n vs. n at the occurring of 
the first maximum of \A\ 2 , at r = 12 A. The other parameters are 8 = and A(0) = 10 -4 . 

FIG. 2. Quantum limit of CARL for p < 1 in the case k = 0. (a) \A\ 2 and (b) (p) vs. r, for p = 0.2, 6 = 5, A(0) = 10" 5 and 
the atoms initially in the state n = 0. We note that (p) = — (2 /p)(|^4_| 2 - |y!(0)| 2 ). 

FIG. 3. Sequential superfluorescent (SF) regime of CARL, (a) \A\ 2 and (b) (p) vs. r, for p = 2, 6 = 0.5, k = 10, and the 
same initial conditions of fig. 2. 
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